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Abstract 

In cosmological scenarios with thermal inflation, extra eras of moduli matter domination, ther¬ 
mal inflation and flaton matter domination exist between primordial inflation and the radiation 
domination of Big Bang nucleosynthesis. During these eras, cosmological perturbations on small 
scales can enter and re-exit the horizon, modifying the power spectrum on those scales. The largest 
modified scale, /cb, touches the horizon size when the expansion changes from deflation to inflation 
at the transition from moduli domination to thermal inflation. We analytically calculate the evolu¬ 
tion of perturbations from moduli domination through thermal inflation and evaluate the curvature 
perturbation on the constant radiation density hypersurface at the end of thermal inflation to de¬ 
termine the late time curvature perturbation. Our resulting transfer function suppresses the power 
spectrum by a factor ~ 50 at /c S> /cb, with k\^ corresponding to anywhere from megaparsec to 
subparsec scales depending on the parameters of thermal inflation. Thus, thermal inflation might 
be constrained or detected by small scale observations such as CMB distortions or 21cm hydrogen 
line observations. 
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I. INTRODUCTION 


Thermal inflation [1-9], a brief low energy inflation that occurs when thermal effects 
hold an unstable flat direction at the origin, occurs naturally in supersymmetric theories 
and is motivated to solve the moduli and gravitino problems [10-12], Additionally, ther¬ 
mal inflation provides a mechanism for baryogenesis [13-20] and has implications for dark 
matter [3, 16, 19], When thermal inflation is included, observable cosmology starts as usual 
with primordial inflation, which generates the density perturbations that go on to form 
galaxies, etc., but then has additional eras of moduli matter domination, thermal inflation 
and flaton matter domination inserted before the usual radiation domination of Big Bang 
nucleosynthesis (BBN). 

During these thermal inflation eras, cosmological perturbations on small scales can enter 
and re-exit the horizon, modifying the power spectrum on those scales, while perturba¬ 
tions on larger scales, corresponding to cosmic microwave background (CMB) or large scale 
structure (LSS) observations, remain outside the horizon, preserving and only modestly 
redshifting their spectrum. 

The observational impact of thermal inflation on the gravitational wave background was 
studied in [21, 22]. Thermal inflation wipes out any potentially observable gravitational 
waves from primordial inflation on solar system or smaller scales [21], but the first order 
phase transition at the end of thermal inflation generates gravitational waves [22, 23] with 
frequencies in the Hz range, though their amplitude may be too small for them to be observed 
in the near future [22]. 

In this paper, we study the impact of thermal inflation on small scale density perturba¬ 
tions with the aim of finding signatures of, or constraints on, thermal inflation. We note that 
a variety of small scale physics, such as ultracompact minihalos or primordial black holes 
[24-26], lensing dispersion of SNIa [27], CMB distortions [28, 30] and the 21cm hydrogen 
line at or prior to the era of reionization [32, 33] could be used as tools for studying the 
effects of thermal inflation on the small scale power spectrum. Several prominent upcoming 
observations are designed for such small scale physics, for example the Primordial Infla¬ 
tion Explorer (PIXIE) [34] and the Polarized Radiation Imaging and Spectroscopy Mission 
(PRISM) [35, 36] for CMB distortions and the Square Kilometre Array (SKA) [37] for the 
21cm hydrogen line. 
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FIG. 1. Characteristic scales of cosmology with thermal inflation. There are four characteristic 
scales, kg,, /cb, kc and kd, corresponding to the comoving scale of the horizon at each of the era 
boundaries, k^, is the largest, and hence most observationally relevant, scale, see Eq. (48). 

In Section II, we outline cosmology with thermal inflation and determine its characteristic 
scales. In Section III, we analytically calculate the thermal inflation transfer function which 
modifies the primordial power spectrum. In Section IV, we summarize our results and briefly 
discuss the possibilities of observing the effects of thermal inflation. In this paper, we set 
h = c = SnG = Mpi = 1. 

II. COSMOLOGY WITH THERMAL INFLATION 

In cosmology with thermal inflation [1, 2, 13-16, 22], there are two different epochs 
of inflation. The first, primordial inflation [38-42], generates the primordial perturbations 
which are observed in the CMB and grow to form galaxies and the LSS. The second, thermal 
inflation, is brief and occurs after primordial inflation but before Big Bang nycleosynthesis, at 
a sufficiently low energy scale to solve the moduli problem, and may affect the perturbations 
on very small scales as we discuss in the following sections. 
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In Fig. 1, we illustrate the thermal inflation eras of moduli domination, thermal inflation 
and flaton domination. They are preceded by primordial inflation plus a possible post¬ 
inflationary era and followed by the radiation domination of BBN. 

For definiteness, we consider a general class of supersymmetry breaking scenarios in which 
supersymmetry is broken in a hidden sector at a scale Mg and transmitted to the observable 
sector via gravitational strength interactions, so that the supersymmetry breaking scale in 
the observable sector is rrig M^/Mpi. We set rrig ~ 10^ GeV. 

In the early universe, the finite energy density, represented by the Hubble parameter H, 
breaks supersymmetry. When H ^ rus, the supersymmetry breaking by the finite energy 
density dominates over the vacuum supersymmetry breaking giving a moduli potential of 
the form 

However, as the Hubble parameter drops through H ~ m$, the moduli potential switches 
to its vacuum form 


with 
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and <Fi — <F 2 ~ Mp\. Thus, at M ~ m$, or t ~ ta in the notation of Fig. 1, the moduli 
start oscillating with Planckian amplitude, dominating the energy density of the universe. 
Furthermore, the moduli have relatively low masses and very weak interactions so this moduli 
oscillation may be sufficiently long-lived to have disastrous effects on, for example, BBN. 
During the era of moduli domination, the moduli abandunce is 


s V 


( 4 ) 


where is the moduli number density and s is the entropy density, but not to spoil BBN 
it should be [43] 

^ < 10-15 to 10-12. (5) 

s 

This is called the moduli problem [10-12]. 

Thermal inflation, which is motivated to solve the moduli problem, is realized by an 
(almost) flat direction with negative mass-squared at the origin, called a flaton. Like the 
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Standard Model Higgs field, the fiaton potential near the origin is 

Unlike the Standard Model Higgs field, but like many scalar field directions in the Minimal 
Supersymmetric Standard Model, the fiaton does not have a stabilizing 0^ term. Instead, 
higher order terms, or the renormalisation group running of the fiaton mass squared, stabilize 
the fiaton potential at a large field value 00 ^ To have zero energy density at the 
minimum, we require Vq ~ and so, for 0o ~ Mpi, we have 

~ \/m^Mpi. ( 7 ) 

Taking 

~ rris ~ 10^ GeV (8) 

gives 

10=^ GeV < Uo^ < 10^^ GeV. (9) 

At the finite temperature of the early universe, the baton’s potential is modified by its 
(assumed) unsuppressed couplings to the thermal bath 

U(0, T) = Uo + ^ - mj) 02 + ..., (10) 

where a is not small. When T > m^/a, corresponding to t ^ tc in Fig. 1, the fiaton is 
trapped at 0 = 0, leaving the potential energy U = Vq which drives the thermal inflation. 
At t ~ tg,, the moduli dominate the universe, but with a comparable amount of radiation 

Pm(ta) ~ Pr(G) > Vq . (11) 

As the universe expands, the moduli and radiation are diluted 

Pm OC a~^ , Pr oc a“^. (12) 

At t ~ fb, the moduli density drops below Vq 

Uo ~ Pm(tb) > Pr(tb) (13) 

and thermal inflation begins. As the universe inflates, the temperature drops 

T oc a~^ (14) 
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and the moduli are diluted to a safely small abundance, see Eq. (24) below. 
At t ~ tc, the temperature drops to the critical temperature 


Tc~ 


a 


(15) 


at which a hrst order phase transition ends thermal inflation [22] and the flaton rapidly rolls 
towards and oscillates about its minimum at 0 = 00 ; leading to a flaton matter dominated 


era 


CX Oj 
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(16) 


The change in state at the end of thermal inflation perturbs the moduli potential by an 
amount 

K) 
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regenerating a moduli abundance 




At 1 ~ td, the flaton decays to radiation at a temperature [16] 


(17) 


(18) 


Td ~ 10^ to 10“^ GeV 


(19) 


and the standard cosmic history of radiation domination, BBN, etc., follows. 

The moduli generated at t ~ ta are diluted by thermal inflation to an abundance [2] 
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where 9*s{T) is the effective number of entropic degrees of freedom at temperature T [44]. 
In models of thermal inflation incorporating a mechanism for baryogenesis [13-16], entropy 
release by Affleck-Dine helds during thermal inflation may typically lead to an effective 
double thermal inflation [15, 16], further diluting the moduli by a factor 
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The moduli regenerated at t ~ tc have abundance [2] 

H in-isl Ti \ Ao^’Oev V/ \ 

s ~ m|Af^, ~ VGeVi t m» ) boiGeV) 


(25) 


Thus for values of in the middle of its expected range, Eq. (9), the moduli can be 
diluted to, and regenerated with, a sufficiently small abundance to satisfy Eq. (5), solving 
the moduli problem. 


A. Characteristic scales 

The thermal inflation eras, illustrated in Fig. 1, determine four characteristic scales, /ca, 
/cb, kc and k^, where 

k^ = (26) 

corresponds to the comoving scale of the horizon at the era boundary We express the k^ 
in terms of the effective thermal inflation parameters 

W, = ln^, (27) 

ax 

which measure the durations of the thermal inflation eras by their number of e-folds of 
expansion, and in turn estimate the N^y in terms of the more fundamental thermal inflation 
parameters Vq; Tc, T^, etc. 
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During thermal inflation, tb < ^ < 4, the energy density p = Vq, so 


kh = e . 


(31) 


At the beginning of thermal inflation, t = th, 
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As discussed above, entropy release during thermal inflation may typically add an extra ~ 6 
e-folds, giving 

A^bc ~ 11. (36) 

More general forms of double or multiple thermal inflation can in principle give even larger 
values of Abe- 


3. kc 


During flaton domination, tc < t < td, the energy density p oc a so 
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At the beginning of radiation domination, t = td, the scale factor is [44] 
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where oq and Tq are the current scale factor and temperature, respectively, and the Hubble 


parameter is 
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therefore Eqs. (26), (40) and (41) give 
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B. The largest characteristic scale 


Eqs. (28), (31) and (37) imply 


kh ^ ■} 


(44) 


/cd kc •) 


(45) 


thus either kh or kd will be the largest scale. Eqs. (31) and (37) give 


kh = e 


d ? 


(46) 


therefore if 

iVbc>^iVed, (47) 

which we will assume, then k^ will be the largest physical scale and hence the one that can 
be most easily observed. Using Eqs. (38), (42) and (46), 
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Modes with k < remain outside the horizon throughout the thermal inflation eras, and 
so are not affected by thermal inflation, while those with k > k^ enter the horizon during 
moduli domination, allowing their evolution to be modifled. Hence, it is expected that there 
could be observable features of thermal inflation at k ^ ky,. 

Modes with k > k^ enter the horizon before moduli domination and so probe that un¬ 
known era, and modes with k > k^ reenter the horizon during flaton domination and so will 
be twice modifled. 

In the next section, we study the evolution of the density perturbations for modes 

ky,< k /Ca, kd . (50) 


III. EVOLUTION OF THE DENSITY PERTURBATIONS 

During the moduli domination and thermal inflation eras, ta < t < tc, we have moduli 
matter (m), thermal radiation (r) and vacuum energy (Vq), with 


P — Pm + Pr + Vq , 
P=^Pr-Vo. 


(51) 

(52) 


To describe the perturbations in moduli and radiation, we define the gauge invariant vari¬ 
ables 
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where Tlsp^ is the curvature perturbation on constant moduli density hypersurfaces and Tlsp^ 
is the curvature perturbation on constant radiation density hypersurfaces. Eq. (A20) gives 
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where q = k/a. 

The physics at t ~ ta is uncertain so we restrict ourselves to ^ t < tc, in which case 


Pv < Pn 


and so Eqs. (55) and (56) reduce to 
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Dehning the curvature perturbation on constant density hypersurfaces and the entropy 
perturbation 

TZsp = ^Usp^ + ^TZsp, ~ Usp^, ( 66 ) 

P P 

Srar = Psp^ ~ Psp-^i (67) 

respectively, Eqs. (58) and (59) are equivalent to 
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where 


- — 2 2 ~ ) , 
Pm + V 3 


which have solution 


n,,(k,t) ^ A„(kA) l + ^H(t)l‘ 'i^^0^ +BUk,t>)§l^ 

S,„r(k,t) = ['P-Si^{k,ti) - 1Zsp,„i^^U)]cos ^ 

+ [n,p,{k,t.)-iispJkA)] 

n-^hnK]aik.t') 


Iti q{k,t') \Jt' 1/3 


G(k,t) = ~g(k,tf 1 A..(*;,ii)ff(i)/^' 


‘ g{k, t') dt' H(t) 


Thermal inflation ends when the temperature of the radiation drops to a critical tem¬ 
perature, given by Eq. (15), triggering a first order phase transition converting the vacuum 
energy Vq into flaton matter. The transition hypersurface matches to constant radiation 
density hypersurfaces before the transition and constant density hypersurfaces after the 
transition. Thus the late time perturbation is given by identifying 

^<5p = TZsp^ {k,t-y (74) 

Using Eq. (62) and taking t\ <tc < U gives 
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The scale 


k\y - (^(tb) H , 


discussed in Section II, is of central importance in this paper, and we precisely define tb, the 
boundary between moduli domination and thermal inflation, by 


a(fb) = 0 
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or equivalently, using Eq. (57), 
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(85) 


14 



where 
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A. Thermal inflation transfer function 

If we take an adiabatic initial condition 


Srar{k, t{) = S^r{k, ti) = 0 (88) 

SO that 

(k, ti) = TZsp, {k, ti) = TZspik, ti ), (89) 

{k, ti) = hsp, {k, ti) = ilsp{k, ti) (90) 


neglect the decaying mode, B^{k,ti) = 0, and take ti ^ 0 and t{ —)■ oo, the effect of the 
thermal inflation era on the curvature perturbation can be expressed as the transfer function 


^ 7^^p(/t, oo) 

nsp{K,0) ’ 


which, from Eq. (85), is given by 


(91) 



and is plotted in Figure 2. On large scales, which remain outside the horizon, the transfer 
function asymptotes to one 
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where, using Eq. (Bl), 
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FIG. 2. The thermal inflation transfer function T{k), where k = k/k\^ and the characteristic scale 
/cb is defined by Eq. (77) and estimated in Eq. (48). The first peak is at (k,T) ^ (1.13,1.21) 
and the first dip is at (3.11,0.269). The analytic form is given in Eq. (92) and the asymptotic 
behaviours are given in Eqs. (93) and (95). 


On small scales, which enter well into the horizon dnring modnli domination and thermal 
inflation, the transfer fnnction oscillates dne to the oscillation of the radiation pertnrbation 
inside the horizon 

T{k) —- cos(i/ifi:) + , (95) 

where we have used Eq. (B4) and 


z/i = 



da 

\Ja{2 + a^) 


2V33y^ 


~ 2.2258. 


(96) 
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At =k/kT^ 


FIG. 3. The thermal inflation transfer function T^^k). 

The power spectrum after thermal inflation is 

P{k) = r\k/h)xP^,,{k), (97) 

where Ppri{k) is the power spectrum of the primordial inflation, see Figure 3. Hence, the 
power spectrum is enhanced by a factor 1.46 at the first peak and, from Eq. (95), suppressed 
by a factor 50 on small scales. 

IV. DISCUSSION 

In this paper, we presented the effect of thermal inflation on the linear evolution of small 
scale density perturbations. By introducing thermal inflation, the post-inflationary history 
is changed to include extra eras of moduli matter domination, thermal inflation and flaton 
matter domination, followed by the usual radiation domination, as in Figure 1. Modes with 
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k > kh enter the horizon during moduli domination before exiting again during thermal 
inflation, and hence are modihed relative to the modes with k < kt, which remain outside 
the horizon. At the end of thermal inflation, the radiation perturbation is converted into the 
hnal adiabatic (curvature) perturbation. The net effect of this for an adiabatic primordial 
perturbation is to suppress the perturbations by a factor ~ 50 on scales /c 3> /cb, see Figures 2 
and 3, and Eq. (95). 

A multicomponent, i.e. adiabatic (curvature) plus entropy (isocurvature), primordial per¬ 
turbation leads to a more complicated, initial condition dependent, result given in Eq. (85). 
Even for modes with k < /cb, which remain outside the horizon, the case with thermal 
inflation can lead to different results compared to without thermal inflation as it is the sub¬ 
dominant radiation perturbation, rather than the dominant moduli matter perturbation, 
that is converted to the final adiabatic perturbation. 

The scale /cb is theoretically estimated in Eq. (48), but the uncertainties are large, with 
anything from megaparsec to subparsec scales being reasonable, and no robust upper or 
lower bound. Current and future observations could provide stronger constraints, or an 
observational signature of thermal inflation. For example, the lack of small scale suppression 
in the primordial power spectrum reconstructed from the recent Planck results [46] suggests 
kb>l Mpc“^ 

The suppression of the power spectrum would reduce the number density of dark matter 
clumps of mass 



where the scale 3/cb is estimated from Figure 3. Depending on the value of /cb, this might 
be able to be seen from 21cm hydrogen line observations such as the Square Kilometre 
Array (SKA), see for example [37], or gamma-ray observations of WIMP-annihilation in, or 
gravitational microlensing of, ultracompact minihalos, see for example [26]. 

Silk damping might be relevant when radiation perturbations enter the horizon during 
either tg, < t < oi t > For tg < t < tc, we estimate the Silk damping scale to be 
much smaller than our characteristic scale k^ [47]. For t > t^, modes near fcb can be affected 
by Silk damping. For 50Mpc“^ ^ /cb ~ lO'^Mpc”^, the dissipation of modes with /c ~ /cb 
would generate a unique pattern of CMB distortions, different from [29-31, 45], due to the 
suppression and oscillations of the thermal inflation transfer function. Future observations 
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such as the Primordial Inflation Explorer (PIXIE) [34] or the Polarized Radiation Imaging 
and Spectroscopy Mission (PRISM) [35, 36] could be used to probe these distortions [45]. 
For 10"^Mpc“^ ^ /cb ~ lO^Mpc”^, the suppression of the power spectrum could diminish 
any effect of the release of energy by Silk damping on Big Bang nucleosynthesis [47]. 
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Appendix A: Review of perturbations in a multicomponent system 


The scalar parts of the perturbed metric and energy momentum tensor are [48] 


ds = (1 + 2A) dt^ — 2B^i dt dx^ — (1 + 277) a^{t) 5ij + 2C 




dx* dx^ 


and 


where ifi = 1/a and 




p = p + 5p , id = -^diV, 


p = p + 5p , = didjTT. 

The Einstein equation gives 

3H (n - ha) +q‘^[n-H(C- 2HC - b) 

il-HA =-H {y + B), 


= \sp, 


d / • 


dt 


(C - 2HC -B)+H(d- 2HC - b) - TZ - A 


= TT, 


7 -1 -1 

-(n- ha) +3H(n- ha) -HA = — 6p+ -q^ir, 

nf \ /V / 9 .“t 


dt V' ■ -- V ' -- -v ---- 2'" ■ 3 

where q = k/a. Taking the derivatives of Eqs. (A4) and (A5) gives 

6p + 3H {6p + 6p) + 2q^H {v + B) - 2H 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 


3n - q‘^ (^C - 2HC - b)] = 0, (A8) 
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7 XI 

- [if (t; + B)] + 3HH {v + B)- -6p + -q\ + if A = 0 

dz ^ Zi o 

corresponding to V ■ T = 0. Decomposing into components 

Sp = 'E.Spx , (p + p)v ^ J2(Px+Px)vx, 


(A9) 




6p = ^5px 

X 


X 

n = '^nx 

X 


(AlO) 


gives 


H {^Px + 3H (ipx + hx) - (px + Px) ivx + B) + (px + px) 37l-q'^(C - 2HC - -S) } = 9. 

(All) 

X] [{Px +px) {vx + B)] + 3H [px +px) (vx + B) + Spx - ^q^'^x + ipx + Px) M ^ 0 

(A12) 

corresponding to J2x V ■ Tx = 0. 

For components which couple only gravitationally, we have V ■ Tx = 0 separately 

Spx {Spx Spx) — q'^ {px Px) {vx B)(pxPx) 3TZ — q^ (C — 2HC — B^ =0, 

(A13) 

d 2 

^ [{Px+Px) {vx + B)] + 3H {px +Px) {vx + B) + 5px--q'^7rx + {px + Px) A = 0. (A14) 

For 


Spx = —Spx , (A15) 

px 

TTx = 0, (A16) 


combining Eqs. (A13) and (A14) and using Eqs. (A4) to (A7) gives 



—Spx)+h(2-3?P-]^^ 


px 

1 + 3^ 

_ px 

p + p+lq^ 


Px J dt 


n 


—dpA+q^— (n-—6px 
Px J Px \ px , 


' .d 

r 


!2sp\ +q^(p + p) 
P 


n 




p ) 

(A17) 


with 


ff 


77- 

'5f^ 

, P 


-6p^Z- 

P X P 


^ if ' 
77 - —6px 
V Px , 


6p" 

P . 


= -E 

X 


q P 


P . 


E if ' 

77 - —6px 
V Px p 


(A18) 

(A19) 
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For two components, X and Y, Eq. (A17) gives 


where 



TZspx ='n — —Spx- 
Px 


(A20) 


(A21) 


Appendix B: Mathematical formulae 

Integrating by parts: 


er^-lf’dpl ^ 

Jo 7 ”^ Jo 


3 

2 pOO 


da 


= 3 

_ 1 
“ 2 


2 + I3^ J J'f ^Ja{2 + q ; 3 ) 

dy /’°° da 


d'y 


0 -yi/2 (^2 + 7^)^^ 7r(27”c^ 7^^^ (2 + 7^) 

f°° da 




n 2 


'0 -^0(2 + a3) 

where we have used 

and 

Integrating by parts: 


[ da(—^ 

lo V2 + q;3 


(]^ry Y 


3/2 


(2 + y 3 ) 3/2 (2 + ^ 3 ) 1/2 

3dy (2 + y=^)^/" 


^5/2 (2 + ^3)1/2 J 


3/2 


foo d-f . 
OK — Sin 
Jo 'y^ 


= — cos 
5 


da 


K 


"! \J a{2 + a^) 




/■°° da 

a 

/■°° da 

•^0 ^a{2 + a^)_ 

•^0 .^ 7(2 + 73 ) 

A a{2 + a^) 


X 




7 + 7 

y3 ; 7o (2 + /33)5/2 


d/3. ' 


= —- cos 
5 


da 




0 ^a(2 + a3)_ 


(Bl) 


(B2) 


(B3) 
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IS f°° , . 

H-/ d'y sm 

K Jo 



da 

^a{2 + a3) 



9 n 

(2 + / 33 ) 5/2 


(B4) 
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